Symmetries in nonlinear Bethe-Heitler process 
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Nonlinear Bethe-Heitler process in a bichromatic laser field is investigated using strong-field QED 
formalism. Symmetry properties of angular distributions of created e~e^ pairs are analyzed. These 
properties are showed to be governed by a behavior of the vector potential characterizing the laser 
field, rather than by the respective electric field component. 
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Quantum electrodynamics (QED) predicts a possibil- 
ity of spontaneous emission of electron-positron pairs 
from a vacuum in the presence of a static electric field []], 
[2I]; the mechanism known as the Sauter-Schwinger mech- 
anism (see, e.g., Ref. [3|). This nonlinear effect was also 
considered for time-dependent electric fields [3, HI • In the 
latter case, the pairs are tunneled through the QED vac- 
uum by a strong low-frequency oscillating field; the point 
is that when an alternating electric field changes slowly 
over the Compton length then it can be treated quasis- 
tationary. In the quasistationary picture, the pairs are 
produced if work performed by the electric field over the 
Compton length provides the energy of 2meC^, with the 
probability determined by the maximum electric field £ 
(above, ttIq is the electron mass). More specifically, the 
probability of pair creation scales as exp(— Tr^cr/^) for 
£ <C £cri which is the Sauter-Schwinger formula with the 
critical electric field ^cr = 1.3 x 10^^ V/cm. 

In this paper, another scenario of pair creation, the so- 
called nonlinear Bethe-Heitler process, in which electron- 
positron pairs are created in collisions of a super-intense 
laser beam with a beam of relativistic targets is inves- 
tigated. We will show that in the case under investi- 
gations, it is the vector potential (not the electric field) 
which determines properties of the process. This is par- 
ticularly interesting in the context of Sauter-Schwinger 
mechanism. In our case, even if we deal with a low- 
frequency laser field, in the reference frame of a target 
particle the field can exhibit, in fact, very rapid oscil- 
lations over the Compton length. Let us remind that 
in usual experiments, incident targets move with very 
high kinetic energies; in the SLAC experiment reaching 
almost 50 GeV that corresponds to the Lorentz factor 
7 10^ [il, 0- Therefore for a typical TiiSapphire laser 
field, the Doppler upshifted frequency of the laser field 
in the target-particle rest frame of reference can reach 
even ttIqC^ . Also the laser field strength in the moving 
reference frame will be significantly enhanced. In partic- 
ular, for record laser fields with intensities of the order of 
(10^^ - 10^^) W/cm^ 0, the corresponding electric fields 
measured in the moving frame of reference can even ex- 
ceed £cv 



In the context of the Sauter-Schwinger mechanism, we 
investigate the role of an incident laser field in the nonlin- 
ear Bethe-Heitler process. This will be done for param- 
eters available today in experiments, that may lead to 
considerable number of e~e+ pairs. Contrary however to 
experiments that have used electrons as target particles 
(see, for instance, in Refs. [il,^), we will use protons. The 
point being that in this case one does not have to account 
for the exchange diagram when calculating the S'-matrix 
amplitude, which makes it more straightforward to in- 
terpret. For our purpose, the angular spectra of product 
particles will be analyzed. The key question is how these 
distributions correlate with symmetry of both the electric 
field and the vector potential characterizing the incident 
laser field. 

We consider the electron-positron pair creation in col- 
lisions of a proton with a bichromatic laser field using the 
theory and numerical methods developed in our most re- 
cent paper [9] . Let us mention that the theory presented 
there, and the current results as well, account for the 
proton recoil. We consider the case when both compo- 
nents of the laser field have commensurate frequencies 
and they are linearly polarized, with the polarization vec- 
tors along the x axis. We also assume that the laser field 
propagates in the opposite z direction. Even though the 
theory introduced in Ref. [9[ is very general, we present 
the results which relate to the reference frame where the 
proton is at rest long before the collision takes place. In 
other words, if p^- and Pe+ stand for the four- momenta 
of created electron and positron, respectively, whereas q\ 
and qi are the incoming and outgoing four-momenta of a 
proton, all numerical results will relate to the reference 
frame in which grf = 0. Consequently, uj will represent 
a Doppler-upshifted fundamental frequency of the laser 
field oscillations. The same approach was used in Ref. [9[. 

To be more specific, we consider here two distinct rep- 
resentations of the four- vector potential which describe a 
two-color laser field. In the first case. 
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A{<j)) = Aosf{4>) , 

/(0) = -^[sin(0)-isin(20) 
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Figure 1. (Color online) In the upper panel, we show shape 
functions for the four- vector potential /((/)) (blue solid line) 
and for the electric field f (red dashed line), describing a 
bichromatic laser field defined by Eqs. {[T]) and {[2]). The laser- 
field parameters are such that cj = meC^ and /i = 1. The lower 
panel presents coarse-grained angular distributions of proba- 
bility rates of created electrons, which have been projected 
onto the xy plane. 



whereas in the second case, 

A{^) = AoeF{^) 
4 



cos((/)) — - cos(2^) 



(3) 
(4) 



with (p = k • X. Here, we have introduced Aq which de- 
scribes the strength of the laser field and is related to the 
relativistically invariant parameter ji = \eAQ\/ {rriQc). As 
anticipated above, the polarization and wave four- vectors 
are e = (0,63^) and k = (cj/c)(l, — e^), respectively. One 
can check that for this configuration, it holds that /c^ = 
and k • A = 0. These conditions constitute foundations 
of the theory developed in Ref. [9|. Moreover, in Eqs. (j2j) 
and dU, we have introduced the so-called shape functions 
/((/)) and F{(j))^ normalized such that 
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and similarly for (F^) = 1/2. As we have discussed in [9|, 
imposing these conditions guarantees that the pondero- 
motive energy of the electron quiver oscillations in a laser 



Figure 2. (Color online) Same as in Fig. [T] but for the bichro- 
matic laser field described by Eqs. Q and (|4|. 



field is the same in both cases. In other words, the pair 
creation threshold energy is also the same, which makes 
it meaningful to compare the rates of pair production for 
both considered cases. Let us note that for our choice of 
the shape functions, it happens that F'{(j)) ~ /(0). This 
means that while the shape function /(</)) describes the 
vector potential for the first choice [Eqs. ^ and (j2j)], it 
relates to the electric component of the laser field for our 
second choice [Eqs. (|3]) and (|4])]. It is the purpose of this 
paper to investigate which of these factors actually deter- 
mines the behavior of probability rates of the nonlinear 
Bethe-Heitler process. 

In the lower panel of Fig. (TJ we present coarse-grained 
angular maps of created electrons projected onto the xy 
plane, for the case described by Eqs. ^ and (j2j) such that 
uj = rrieC^ and /i = 1. These results have been obtained 
using the Monte-Carlo method that we developed in [lo|. 
Here, 3 x 10^ sample points were collected, which then 
were averaged according to the prescription described in 
Ref. [lo|- We have also calculated the total probability 
rate of pair creation which turned out to be 0.11 x 10~^ in 
relativistic units, with an estimated error less than 1%. 
It is crucial to note that in the considered case, electrons 
(positrons as well) are created symmetrically with respect 
to the X axis, i.e., with respect to the polarization vector 
s. In the upper panel of Fig.[TJ we demonstrate the shape 
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Figure 3. (Color online) Angular maps of i?4(gf,Pe-) (1^^ 
column) and i?7(gf , p^- ) (right column) in the case when Qf = 
—2mecez. The parameters of the laser field have been chosen 
such that cj = rrieC^ and /i = 1. While the upper panels 
refer to the vector potential specified by Eqs. ([1]) and ([2|), the 
lower panels are for the vector potential described by Eqs. ^ 
and dH). For visual purposes, the rates have been raised to 
power 1/5. 



Figure 4. (Color online) Angular maps of Ri83{Qi,Pe-) (left 
column) and RiseiQi^Pe-) (right column) in the case when 
Qi = —2mecez. The parameters of the laser field are uj = rrieC^ 
and /i = 10. The upper panels correspond to the case specified 
by Eqs. ([T]) and (|2]), whereas the lower panels are for the vector 
potential described by Eqs. ^ and (|4]). For visual purposes, 
the rates have been raised to power 1/10. 



functions f{(j)) and /'(0) which describe either the vector 
potential (blue solid line) or the electric field (red dashed 
line), respectively. There is a clear symmetry of f{(j)) as 
opposed to an asymmetry of /'(^). One can see therefore 
that the symmetry of the vector potential is reflected in 
the pattern of angular distributions of created particles. 

Fig. [2] relates to the second choice of the vector po- 
tential [see, Eqs. (j3j) and (j4])], with the same laser-fleld 
parameters as in Fig.[TJ In this case, the total probability 
rate of pair creation is 0.46 x 10~^ in relativistic units, 
with an error less than 1%. In the upper panel, we plot 
the shape functions F{(j)) and f {(/)). As we have already 
mentioned, in the current case F{(j)) describes the vector 
potential (blue solid line) whereas f{(j)) corresponds to 
the electric fleld (red dashed line). This time, the vector 
potential shape function F{(j)) is asymmetric while the 
electric fleld shape function f{(j)) stays symmetric. In 
the lower panel of Fig. O we present the coarse-grained 
angular distributions of produced electrons. Clearly, an 
asymmetric pattern is observed there. Again, the behav- 
ior of the electric fleld component seems to be irrelevant 
when it comes to properties of pair creation. The same 
will be elaborated below, when analyzing partial angular 
distributions of created particles. 

In the following, R^^^ {Qi^Pe- ) denotes the partial triply 
differential probability rate of pair creation when the 
value of energy transfer from the colliding proton as well 
as orientations of the proton and the emitted electron 
are flxed. In Refs. [ol, lU, |T2|, we have also used the 

quantity RN{Qi^Pe-) = "^Zi ^N^i^^^Pe-) which is a sum 



of the partial differential probability rates for the same 
four- momenta conservation condition [III, [HI- Thus in 
Figs. [3] and m we present angular maps of RN{Qi,Pe-) 
in function of the electron polar and azimuthal angles, 
Oq- and (Pq- , for chosen numbers of photons N absorbed 
from the laser fleld (for more details, see Ref. @). 

In Fig. [31 we demonstrate angular maps of the triply 
differential probability rates of pair production by 4- 
laser photons (A^ = 4), R^iqi^Pe-)^ and by 7-laser pho- 
tons (TV = 7), Rjiqf^Pe-)' The results are for the case 
when the Doppler up-shifted frequency of the fundamen- 
tal color fleld is = mgC^, = 1, and the proton recoil 
is grf = —2meCez. While the upper panels of Fig. [3] are 
for the shape function /(</>) [Eq. ([2])], the lower panels 
are for the shape function F{(j)) [Eq. (|4j)]. We observe 
that these angular maps are stretched toward polar an- 
gles 6q- close to TT, which is in the propagation direction 
of the laser fleld. This has been realized to be a gen- 
eral feature of strong-fleld multiphoton processes in the 
relativistic domain, where the motion of electrons in the 
direction of the laser fleld propagation is driven by the 
magnetic component of the Lorentz force MM- We 
note also that these angular maps exhibit the same type 
of behavior as the vector potential does, namely, they 
are symmetric when the vector potential is symmetric 
(upper panels), and asymmetric otherwise (lower pan- 
els). In contrast to these very few photon processes, in 
Fig. m we present maps for RissiQf^Pe-) = 183) and 
RiseiQf^Pe-) (A^ = 186), as denoted below each panel. 
This time, the laser fleld is even stronger than in Fig. [3l 
since /i = 10. Even though in the current case, we ba- 
sically reach the quasistationary regime of pair creation. 
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that is believed to occur when /i > 1 [Tol, we 
stih observe the same symmetry properties of angular 
maps of created electrons as in Fig. [H These proper- 
ties are governed by the behavior of the vector potential, 
not bythe electric field. If the Sauter-Schwinger mecha- 
nism [l|,[2| was responsible for the nonlinear Bethe-Heitler 
process this would not be the case. At this point, let us 
also mention that for this strong field (corresponding to 
/i = 10) we reach the limit of computational capacity in 
our Monte Carlo simulations. For this reason, we do not 
present the respective result for the total probability rate 
of pair creation. 

Let us note that angular distributions of created parti- 
cles through nonlinear Bethe-Heitler scenario have been 
studied before for the case when an incident laser field 
is treated as a monochromatic plane wave (see, for in- 
stance, Refs. [nl, [13 )• In this case, the calculated an- 
gular spectra turn out to be symmetric under a mirror 
reflection with respect to the polarization direction of the 
laser field. At the same time, both the vector potential 
and the electric field are symmetric in the plane wave 
approximation. By adding the second color field, one 
can break the symmetry of either the vector potential or 
the electric field, and analyze how this may infiuence the 



pair production. This is the essence of the present paper 
and also a suggestion for experimentalists to validate our 
idea. 

In summary, symmetry properties of nonlinear Bethe- 
Heitler process by a two-color laser field have been inves- 
tigated. These properties have been proved to depend on 
the vector potential characterizing the laser field, not on 
the electric field component. This can be understood if 
we note that in the strong-field QED formalism which is 
based on the so-called Volkov waves [ll|, [l8| , the maxi- 
mum of probability rates of pair creation is achieved for 
saddle-point solutions of the classical action describing a 
charged particle in the laser field. In our case, these sta- 
tionary points are approached when produced particles 
are detected with asymptotic momenta such that their 
component along the polarization direction is equal to 
eA. Hence, the vector potential does determine proper- 
ties of ultrastrong pair creation. In this sense also, the 
Sauter-Schwinger mechanism which relies on the electric- 
field dependence of pair creation probability rates ap- 
pears to be irresponsible for the process under investiga- 
tions, at least for the chosen laser-field parameters. 
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